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Fluctuations in the two-dimensional one-component plasma and 
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We study the distribution of the mean radial displacement of charges of a 2D one-component plasma in the 
thermodynamic limit A —>■ oo at finite temperature /3 > 0. We compute explicitly the large deviation functions 
showing the emergence of a fourth-order phase transition as a consequence of a change of topology in the plasma 
distribution. This weak phase transition occurs exactly at the ground state of the plasma. These results have 
been compared with the integrable case (finite N) of plasma parameter = 2. In this case the problem can be 
mapped to the stationary properties of 2D Dyson Brownian particles and to a non-Hermitian matrix model. 


Introduction. The two-dimensional one-component plasma 
(2D-OCP) is one of the most basic statistical mechanics fluid 
(i.e. non lattice) models in dimension greater than one. Very 
recently, Sandier and Serfaty [1] have rigorously established 
a connection between the Ginzburg-Landau model [2] for su¬ 
perconductors in the critical regime with vortices and the 2D- 
OCP. In this regime the interaction between vortices is well- 
described by a one-component plasma in a quadratic confining 
potential. Heuristic arguments suggest that the very same ap¬ 
proximation could be valid to describe more general vortex 
systems, most notably in superfluids or Bose-Einstein conden¬ 
sates [3, 4]. 

The 2D-OCP (also known as log-gas) is a system of clas¬ 
sical pointlike particles of same charge q (one species of par¬ 
ticles) immersed in a two-dimensional domain with a neu¬ 
tralising background of opposite charge. The interaction be¬ 
tween the charges is (1/2) — r)|), where the 

Coulomb potential is solution of the 2D Poisson equation, 
v{f) = — log(|f^/L). The Boltzmann-Gibbs canonical mea¬ 
sure of a 2D-OCP of N particles at inverse temperature /3 in a 
quadratic potential is 

P/3 ({4}) = (1) 

k 

where = {xk, j/fc) G is the position of the fc-th particle 
(fc = 1,..., N), Tk = |rl| and rij = jf) — fj\. The thermody¬ 
namic state of the system is characterised essentially by the 
two-dimensional plasma parameter 7 = Hereafter the 
length L that Axes the zero of the potential and the elementary 
charge q will be scaled to unity. (The plasma parameter 7 is 
therefore identified with the inverse temperature /3.) 

Besides being an approximation of the Ginzburg-Landau 
model in the vortices phase, these log-gases are interesting per 
se, being ubiquitous in many fields of physics and mathematics. 
At the special value /3 = 2 of the plasma parameter, the model 
is exactly integrable and deserves particular attention. In this 
case, (1) is the stationary distribution of a 2D Dyson Brown¬ 
ian motion [5, 6], namely a Brownian motion of N particles 
with logarithmic repulsion in the plane confined in a harmonic 


potential. Again at /3 = 2, with the identification ~ C, (1) 
is the joint distribution of the eigenvalues Xk = Xk + ipk of 
non-Hermitian random matrices with iid complex Gaussian 
entries, the Ginibre ensemble [7]. There has been extensive 
work on 2D-OCP systems in the mean-field theory; numerical 
studies are also available, while exact results can be derived in 
the integrable case /? = 2 [8-17]. 

In this paper we present a study on the large deviations 
(rare events) of 2D-OCP systems. We will consider a rather 
simple observable, namely the mean radial displacement 
An = N~^ Tk of the charges of the plasma (the vor¬ 
tices of the Ginzburg-Landau model in the critical phase). 
Here we compute explicitly the large N limit probability law 
= {S{x — An)) of An and its rescaled Laplace 
transform (see equations (3)-(5) 

and (4)-(7) below) for all /? > 0 (the angle brackets denote 
the canonical average). As a byproduct, at /3 = 2 these results 
provide large deviations functions for the Ginibre ensemble 
and 2D Dyson Brownian particles. 

Results on the probability of rare events are popular in the 
one-dimensional case, in particular in the literature on Hermi- 
tian random matrices [18-25]. In most of the models analysed 
the 2D-OCP is confined on the line and undergoes a phase 
transition, usually of third-order if the plasma is constrained 
by a hard wall (for an interesting review see [26]). Large de¬ 
viations have been much less explored in two dimensions and 
few results are available [27-30]. Here we show explicitly that 
the statistical distribution of the mean displacement A^v in the 
plane unveils a fourth-order phase transition of the 2D-OCP re¬ 
lated to a change of topology in the equilibrium distribution of 
the plasma. This phenomenon is absent in one dimension. We 
point out that the extraordinary weak singularity of yN,p{x) 
(or yN,p{s)) occurs at the average value a: = (A^r) (at s = 0 
in the Laplace domain), i.e. at the ground state of the 2D 
plasma. Therefore this phase transition may be detected by the 
anomalous behaviour of the high order cumulants of Ajy. 

Results. The general theory of large deviations of log- 
gases [29, 30] ensures that the Laplace transform CPAr^^(s) of 
A N has a thermodynamic limit 

J(s) = - \im^ ^ log Tn.p {s) (3) 
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FIG. 1: (color online) Top: Rate function (['(a;) from (20) and its 
leading terms for x —>■ 0, x ~ 2/3 and x —> oo (see (7)). At 
X = (Aiv) = 2/3 the function is not analytic. Central: Support E(s) 
of the equilibrium distribution /rj of the plasma for different values 
of s. From left to right s = 1, 0, —0.5. Bottom: The corresponding 
radial distributions d/r/(|r|, ip). 


for all s S K. Such a limit is the scaled cumulant generat¬ 
ing function (cumulant GF) of A at. Then, by the Gartner- 
Ellis theorem [31, 32], the probability law of Ajy satis¬ 
fies a large deviation principle [48] with speed i.e. 

yN,p{x) ~ 'i'(a:)^ where the rate function T'(a:) is the 

Legendre-Fenchel transform of the cumulant GF (3): 

^{x) = - lim log Tat,;? (a;) = inf [J(s) - sx]. (4) 

Af->oo pl\^ s 


The challenge here is to compute explicitly the large devia¬ 
tion functions J(s) and 'I'(x). We will show that the cumulant 
GF has the following expression independent of /3: 


J(s) = Jarsinhf^)-^ + ^ 

^ ^ 2 V2/ 4 48 

(s^-FlO) \/s2-F4-|s|3 


(5) 


This is the central result of the paper. Before turning to the 
derivation, we make few observations and we discuss the con¬ 
sequences of (5). The cumulant GF is strictly concave and 
satisfies the normalization condition J(0) = 0. Notice that 
the large deviation function is not analytic. Tht fourth deriva¬ 
tive of J(s) is discontinuous at s = 0 with a finite jump 
j(4)(o-) _ j(4)[q+] = X; this is a consequence of difourth- 
order phase transition in the equilibrium distribution of the 


2D-OCP. Phase transitions of 2D-OCP constrained on a line, 
i.e. of eigenvalues of Hermitian random matrices, are well 
understood. In two dimensions, however, they are more gentle 
than for Hermitian matrices, where the singularities are of third 
order. This behaviour corresponds to the physical intuition that 
the logarithmic repulsion has stronger effects in ID than in 
2D gases. Here the singularity emerges in the thermodynamic 
limit as a consequence of a change in topology (disk-annulus) 
of the support of the equilibrium distribution of the charges. 
This mechanism is absent in one dimension. We note that 
the phase transition would not occur in the absence the log¬ 
arithmic interaction, since without it (2) is the Hamiltonian 
of a system of free particles, which cannot undergo a phase 
transition — this case is discussed in more detail in the last 
section. As a final remark, we stress that the phase transition 
occurs at the ground state (s = 0) and hence is not related 
to the atypical fluctuations (this is another departure from the 
one-dimensional scenario for continuous statistics). 

At leading order in N, the hrst three cumulants of Aat can 
be obtained by differentiation of J(s) at s = 0 

2 1 1 
«i(AAr) = K2{^n) = K3(Aiv) = (6) 


At /? = 2 the hrst two cumulants (average and variance) agree 
with asymptotic formulas for the Ginibre matrix ensemble [33]. 
The nonzero third cumulant measures the deviation from the 
Gaussian behaviour at leading order in N. Higher cumulants 
cannot be obtained by differentiation of J(s) (we will elaborate 
more on this point later). Since J(s) is everywhere differen¬ 
tiable, we can apply the Gartner-Ellis theorem and recover 
the rate function 'I'(x) using (4). The limiting behaviour of 
yN.p{x) « may be summarized as follows 


J’Ar,/3(x) 


exp 

1(f) 


-(2/3)/3Ar^(a:"-4/9) 

for X <C 1 



for X ~ 1 

-(l/2)/3Af^(a:;"-4/9) 

for X 1 


(7) 


The rate function 41 (x) is plotted in the top panel of Eig. 1. 
Erom the Gaussian form around x ~ 2/3, one easily reads 
off the values for mean and variance: (Aat) = 2/3 and 
var(AAr) = l/(2/3A^), according to (6). On the other hand, 
the rate function is not quadratic and this implies the polyno¬ 
mial and sub-Gaussian far tails of yN,p(x) for x —?► 0 and 
X —>■ oo, respectively, in (7). 

We stress that the exact computation of J(s) is essential 
to unveil the phase transition. Standard computations usually 
concern with the quadratic approximation about the ground 
state of the 2D-OCP (Gaussian approximation of J{s) about 
s = 0). However the fourth-order singularity at s = 0 does not 
emerge at the level of quadratic (i.e. 0(s^)) effects. 

Derivation. In order to compute the cumulant GE (3), 
we cast the Laplace transform of Aat as a ratio of two 
partition functions T’a/./ 3 (s) = where Zn,p{s) = 
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J g jjjg partition function of the 2D-OCP 

with modihed energy 


the inner and outer radii of the support S(s). From (14), the 
equilibrium distribution of the charge has the explicit form 




Vs{r) 


+ N'^Vsirk) ( 8 ) 

k 

1 9 

-r‘^ + sr . (9) 


dM*(rO 


+s \ 
V 27r|f1 ) 


If ei:(s) dr, 


ro(s) =max{0,-s}, Ro{s) 





(15) 

(16) 


Clearly, at zero temperature /3 —oo the only contribution 
to Zn, 0 {s) comes the configurations of the system that min¬ 
imize the energy iF(fi,..., fV; s). At finite temperature 
/3 < oo in the large N limit, the picture is not very dif¬ 
ferent in the sense that the‘ typical’ configurations of the 
system are still the minimizers of the energy and one can 
characterize the probability of ‘atypical’ configurations. For 
large N, we have H{ri,... ,f)v; s) = iV^(iFs[/iAr] -f o(l)) 
where fiN = Ihe normalized distribution of the 

charges and Ffs[-] is the mean-held energy functional: 

Hs[lA = -\ jj A^J.{f)d^x{?)\og\r-?\+ y d/x(f)14(|f1). 

( 10 ) 

These considerations justify a saddle-point evaluation of the 
large deviation function J(s) as excess free energy 




J g-/3ir(fi,...,fjv;s) 
J g-/Sfr(fi,...,fiv;0) 


( 11 ) 


where /r* (resp. /rj) is the minimizer of the energy functional 
Hs[-] (resp. iTo[-]): 


iFs[Ms] = M > 0, y ^ d/r(fO = l| • (12) 


Hence, the problem of computing J(s) reduces to hnd the equi¬ 
librium conhguration n* of the charges of plasma in the large 
N limit. Find the conhguration of the 2D-OCP that minimizes 
a mean-held energy functional like (10) is one of the central 
objectives in potential theory [34]. Fortunately, minimization 
problems such as (12) can be solved in closed form with a 
judicious use of classical electrostatics. We hrst observe that 
the problem is radially symmetric and the radial potential Vs is 
convex. Therefore, the equilibrium measure inherits the radial 
symmetry of the energy d/x*(f) = (27r)“^d/i*(r)dv5 and is 
supported in general on the annulus 

i:(s) = {re ro(s) < r < Ro{s)} , (13) 

where we have switched to polar coordinates r = (r, (p). At 
equilibrium, the Gauss law on a generic domain of the plane 
holds. For a centered disk of radius r, the Gauss law reads 

V'{r)r = 2Tr j d/i*(r')r'. (14) 

Jo 

Using the positivity and the normalization of p,*, it is easy 
to hnd the conditions Vj,'(ro) = 0 and Vs{Ro)Ro = 1 on 


We see that the optimal distribution of the plasma experiences 
a change of topology driven by s: p* is supported on a disk 
for s > 0 and on an annulus for s < 0 (see Fig. 1). The fourth- 
order singularity of J(s) at s = 0 witnesses this sudden disk-to- 
annulus change of topology. In particular, for s = 0 we get the 
uniform measure on the unit disk: d/iQ(f^ = 2 l|f|<idr. This 
corresponds to the circular law [7, 35-37], one of the most 
celebrated results on the empirical density of the eigenvalues 
of non-Hermitian random matrices. For a generic value s € K, 
the distribution p* is the typical conhguration of the 2D-OCP 
with a hxed value of A^r given by 


a;(s) 


dM:(f)|r1 


(s^ -I- 4)^^^ — 6s — |sp 

12 


(17) 


The energy of the equilibrium conhguration at s = 0 (i.e. at the 
ground state) is Hq[pq] = 3/8, a value known as id-stability 
bound of the 2D-OCP [38, 39]. For a generic value of s, the 
minimal energy attained at p* is 


HsIp:] = 


+ VsiRo) - logi?o 


(18) 


Evaluating (18), we get the cumulant GF J(s) in (5) as excess 
energy (11). 

Since J(s) is strictly concave and differentiable, using the 
properties of the Legendre-Fenchel transform, we can write 
the thermodynamical relation J(s) — 'I'(a;) = sx, where the 
conjugate variables x and s are related by 


x{s) = J'{s) and s{x) = —T''(x), (19) 


with x{s) given in (17). These relations, supplemented with 
the conditions J(0) = T'(a;(0)) = 0, provide the following 
simple expression for the rate function T'(a;): 

'^{x) = — f s{x')dx', (20) 

•'I 

where s{x) is the (real non increasing) inverse of (17) and we 
have used a;(0) = 2/3. One may verify that the cumulant 
GF (5) can be computed also from (19). The asymptotic analy¬ 
sis of T* (x) in (20) provides the limiting behaviour summarized 
in (7). 

The integrable case /? = 2. As already mentioned, the value 
/3 = 2 is special. At this value of the plasma parameter the 
2D-OCP canonical measure has a determinantal structure and 
the partition function is integrable at finite N. This fact has 
been largely exploited in random matrix theory. For /3 = 2, 
the Laplace transform of radially symmetric linear statistics 
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FIG. 2: (color online) Top: The cumulant generating function J{s) 
in (5) (solid line) compared with the exact finite N formula (21) 
(at P = 2) for = 16 (dots). Already for such a small value 
of N, the large deviation function J{s) matches the finite N ex¬ 
pression of —1/(2^"^) log Tiv,^^2(-5) with an error 0(10“^) in the 
plotted region. Bottom: Rescaled cumulants for 

m = 1, 2,3,4 at increasing values of N for (3 = 2. The first three 
cumulants have a regular behaviour in N and approach their leading 
order (6), while N^K 4 behaves irregularly and does not converge for 
large N. 


has been computed exactly by Forrester [33]. By specializing 
his result to we obtain the following finite N formula at 
(crucially) ^ = 2: 


N 


5PAr,/3=2(s) = 11 


rdte-V-i 


( 21 ) 


In Fig. 2 we show that our large N result (5) agrees with this 
finite N result even for moderate values of N. We stress how¬ 
ever that the large deviation formula is valid for any value of 
/3 > 0. Moreover, the route presented here shows explicitly 
that, even at /3 = 2, the function — l/(/37V^) log J’Ar_/ 3 (s) an¬ 
alytic at finite N (see (21)) develops a non-analyticity as a 
result of the thermodynamic limit. To understand better the 
emergence of the singularity we have performed another nu¬ 
merical experiment. For (3 = 2 the (unordered!) collection 
of displacements {r^} of the charges is distributed as a col¬ 
lection of independent random variables where 

> 0 has density fk{x) = (2/r(fc))a;^^“^ exp(—a;^) and 
moments = r(fc -|- m/2)/T{k) [40, 41]. The cumulants 
Km ijn > 1) are linear functionals for independent variables 
and hence we have Kmi^N) = N~^ '^k=i i^m{^k)- We 
have evaluated numerically this sum for increasing values of 
N and compared the output with our N ^ 1 results on the 
cumulants (6) (see Fig. 2). For m = 1, 2,3, the sum asymptoti¬ 
cally approaches (6), while the fourth cumulant K 4 ^{Apf) shows 
a bizarre behaviour for large N. For (3 = 2, this is a further 


confirmation of the weak singularity of the large deviations 
functions. 

Role of the logarithmic interaction. In order to convey the 
key ideas of this paper, we conclude with a brief digression to 
clarify the role of the 2D Coulomb interaction in the mecha¬ 
nism of the phase transition. The computation of J(s) led to 
study the minimisation problem (12). For s > 0 the effective 
external potential 14 (r) in (9) is convex and the distribution 
of charges p* is supported on a disk. When s is negative the 
convexity of 14 (r) is broken and the plasma distribution con¬ 
centrates on an annulus. One might think that this change in 
the potential is sufficient to modify the topology of the charge 
distribution and to induce a phase transition even in absence of 
interaction between the charges. This is incorrect, as we will 
demonstrate. 

Removing the logarithmic interaction from (2) turns the 
plasma into a non-interacting gas with single-particle energy 
H ({Dc}) = (1/2) inverse temperature (3 > Q [49]. 

Since the particles are independent, the single-particle distribu¬ 
tion is Gaussian (the system is no more confined in a bounded 
region) and A^v = J^k i® ^ independent and 

identically distributed displacements. Hence, an elementary 
computation shows that (A^v) = •\/7r/(2/3), var(A 7 v) = 
(4 — tt)/{2(3N) and that the Laplace transform is Tjv,;S(s) = 
(g-/3s7VAN) = [1 — ^y2^^P se^p{—(3s'^/2)^{—l3s)]^ , where 
$(a;) is the standard Gaussian cumulative function. (Note 
the different scaling in N and /3-dependence in the above ex¬ 
pressions compared with the 2D-OCR) Therefore, in absence 
of interaction the free energy of the system is analytic in s: 
non-interacting particles do not undergo phase transitions. 

Conclusions. Many of the ideas presented here can be gener¬ 
alized to other problems. In particular, the method to compute 
large deviations formulas for 2D-OCP is applicable to a large 
class of radially symmetric observables Am = /(l^ I) "'ith 

/(r) continuous. The problem of discontinuous statistics, like 
the charge fluctuations (i.e. the index problem in random ma¬ 
trices), is delicate and needs more refined tools (see [42] for 
recent developments). 

We detect a novel and unusual type of weak phase transi¬ 
tion due to a change of topology in the plasma distribution. 
Therefore this work naturally prompts the following question: 
is this fourth-order phase transition on the plane as ubiquitous 
and universal as the third-order transition for charged particles 
on the line? We mention that similar disk-annulus phase tran¬ 
sitions have been detected for non-Hermitian non-Gaussian 
random matrices [43, 44]. However, non-Gaussian ensembles 
do not have a 2D-OCP interpretation and it is not clear how to 
relate these works with our findings. 

An interesting question is to look at the finite N corrections 
of the large deviations principle (3)-(4). At microscopic scales 
it has been shown [1, 17, 39] that, among all lattice configura¬ 
tions, the charges of the plasma prefer to arrange themselves 
in a triangular one (the very same of the Abrikosov lattice [45] 
of vortices in the Ginzburg-Landau model). Showing that the 
triangular lattice is the global minimizer of the 2D-OCP energy 
requires a careful analysis beyond the mean-field limit. Such a 
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conjecture is still unsolved even in the integrable case /3 = 2. 
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